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AN EVALUATION APPROACH TO COMPUTING 
INVARIANTS RINGS OF PERMUTATION GROUPS 



3 NICOLAS BORIE AND NICOLAS M. THIERY 

o 

■^-J Abstract. Using evaluation at appropriately chosen points, we propose a 

Grobner basis free approach for calculating the secondary invariants of a finite 
C'~'3 permutation group. This approach allows for exploiting the symmetries to 

confine the calculations into a smaller quotient space, which gives a tighter 
control on the algorithmic complexity, especially for large groups. This is 
confirmed by extensive benchmarks using a Sage implementation. 



1. Introduction 



o 

5 Invariant theory has been a rich and central area of algebra ever since the eigh- 

1^ 6 teenth theory, with practical applications [DK02[ § 5] in the resolution of polyno- 

C3 r mial systems with symmetries (see e.g. [Col97aj . |Gat90| . |Stu93[ § 2.6], |FR09j ). 

^ 8 in effective Galois theory (see e.g. |Col97b) . |AbdOO| . IGKOO) ). or in discrete math- 

i— ' 9 ematics (see e.g. [ThiOOl IPTOI] for the original motivation of the second author). 

^ 10 The literature contains deep and explicit results for special classes of groups, like 

^ 11 complex reflection groups or the classical reductive groups, as well as general re- 

12 suits applicable to any group. Given the level of generality, one cannot hope for 

13 such results to be simultaneously explicit and tight in general. Thus the sub- 
OO 14 jcct was effective early on: given a group, one wants to calculate the properties 

15 of its invariant ring. Under the impulsion of modern computer algebra, computa- 

16 tional methods, and their implementations, have largely expanded in the last twenty 
^ 17 years |Kem93[ IStu93[ IThiOll IDK021 IKin07b[ IKinOTa) . However much progress is 

18 still needed to go beyond toy examples and enlarge the spectrum of applications. 
• . 19 An important obstruction is that the algorithms depend largely on efficient com- 

^ ^ 20 putations in certain quotients of the invariant ring; this is usually carried out using 

^ 21 elimination techniques (Grobner or SAGBI-Gr5bner bases), but those do not be- 

22 have well with respect to symmetries. An emerging trend is the alternative use 

^ 23 of evaluation techniques, for example to rewrite invariants in terms of an existing 

24 generating set of the invariant ring j GST06[IP5W09] . 

25 In this paper, and as a test bed, we focus on the problem of computing 

26 secondary invariants of finite permutation groups in the non modular 

27 case, using evaluation techniques. 

28 In Section [2] we review some relevant aspects of computational invariant the- 

29 ory, and in particular discuss the current limitations due to quotient computations. 

30 In Section |3j we give a new theoretical characterization of secondary invariants in 

31 term of their evaluations on as many appropriately chosen points; this is achieved 

32 by perturbating slightly the quotient, and using the grading to transfer back results. 

33 In Section |4j we derive an algorithm for computing secondary invariants of permu- 

34 tation groups. We establish in Section [5] a worst case complexity bound for this 
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algorithm. This bound suggests that, for a large enough group G, at least a factor 
of \G\ is gained. This comparison remains however sloppy since, to the best of our 
knowledge and due to the usual lack of fine control on the complexity of Grobner 
bases methods, no meaningful bound exists in the literature for the elimination 
based algorithms. Therefore, in Section [6] we complement this theoretical analysis 
with extensive benchmarks comparing in particular our implementation in Sage 
and the elimination-based implementation in Singular's .GPS98,, Kin07bJ . Those 
benchmarks suggest a practical complexity which, for large enough groups, is cubic 
in the size ?t.!/|G| of the output. And indeed, if the evaluation-based implementa- 
tion can be order of magnitudes slower for some small groups, it treats predictably 
large groups which are completely out of reach for the elimination-based imple- 
mentation. This includes an example with n = 14, \G\ — 50,803,200, and 1716 
secondary invariants. 

We conclude, in Section[7j with a discussion of avenues for further improvements. 

2. Preliminaries 

We refer to [5ta79l IStu93l ICLU971 ISmi971 IKem98l IDK02j for classical literature 
on invariant theory of finite groups. Parts of what follows are strongly inspired 
by [Kem98| . Let V be a K- vector space of finite dimension n, and G be a finite 
subgroup of GL{V). Tacitly, we interpret G as a group of rt x n matrices or as a 
representation on V . Two vectors v and w are isomorphic, or in the same G-orbit 
(for short orbit), li a ■ v — w for some a d G. 

Let X := {xi, . . . ,a:„) be a basis of the dual of V, and let K[x] be the ring of 
polynomials over V. The action of G on y extends naturally to an action of G on 
K[x] by (T • p :~ p o a~^. An invariant polynomial, or invariant, is a polynomial 
p e K[xi, . . . , Xn] such that a - p = p for all a G G. The invariant ring K[x]'^ is the 
set of all invariants. Since the action of G preserves the degree of polynomials, it is 
a graded connected commutative algebra: K[x]'-^ ~ 0^>pK[x]^, with K[x]q^ « K. 
We write K[x]^ = 0^^pK[x]^ for the positive part of the invariant ring. The 
Hilbert series of IK[x]'^ is the generating series of its dimensions: 

oo 

H{K[xf,z) :=^z'^dimK[x]^. 

d=0 

It can be calculated using Molien's formula: 

This formula reduces to Polya enumeration for permutation groups. Furthermore, 
the summation can be taken instead over conjugacy classes of G, which is relatively 
cheap in practice. 

A crucial device is the Reynolds operator: 

R: K[x] — > K[:>if 

which is both a graded projection onto K[x]'-' and a morphism of K[x] '^-module. 
Note that its definition requires charK not to divide |G|, which we assume from 
now on (non-modular case). 
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72 Hubert's fundamental theorem of invariant theory states that K[x]'^ is finitely 

73 generated: there exists a finite set S of invariants such that any invariant can be 

74 expressed as a polynomial combination of invariants in S. We call S a generating set. 

75 If no proper subset of S is generating, 5 is a minimal generating set. Since K[x]'^ 

76 is finitely generated, there exists a degree bound d such that Kfx]*^ is generated by 

77 the set of all invariants of degree at most d. We denote by /3(K[x]'^) the smallest 

78 degree bound. Noether proved that /3(K[x]'^) < \G\. 

79 Thanks to the grading, for AI a set of homogeneous invariants, the following 

80 properties are equivalent: 

81 (i) M is a minimal generating set for K[x]'^; 

82 (ii) M is a basis of the quotient K[x]'^/K[x]^ . 

83 Therefore, even though the generators in M are non canonical, the number of 

84 generators of a given degree d in M is: it is given by the dimension of the component 

85 of that degree in the graded quotient K[x]'^/K[x]'j^ . There is no known algorithm 

86 to compute those dimensions, or even just /3(K[x]'^), without computing explicitly 

87 a minimal generating set. 

88 The previous properties give immediately a naive algorithm for computing an 

89 homogeneous minimal generating set, calculating degree by degree in the finite 

90 dimensional quotient up to Noether's bound. There are however two practical 

91 issues. The first one is that Noether's bound is tight only for cyclic groups; in 

92 general it is very dull, possibly by orders of magnitude. The second issue is how to 

93 compute efficiently in the given quotient. We will get back to it. 

94 By a celebrated result of Shepard, Todd, Chevalley, and Serre, K[x]'^ is a poly- 

95 nomial algebra if and only if G is a complex reflection group. In all other cases, 

96 there are non trivial relations (also called syzygies) between the generators; how- 

97 ever K[x]'^ remains Cohen-Macaulay. Namely, a set of m homogeneous invariants 

98 (01, ... , 9n) of K[x]'^ is called a homogeneous system of parameters or, for short, a 

99 system of parameters if the invariant ring K[x]'^ is finitely generated over its subring 

100 K[0i, . . . , On]. That is, if there exist a finite number of invariants (ryi, . . . , rjt) such 

101 that the invariant ring is the sum of the subspaces rii.K[9i, . . . ,6'„]. By Noether's 

102 normalization lemma, there always exists a system of parameters for K[x]'^. More- 

103 over, K[x]'^ is Cohen-Macaulay, which means that K[x]"-^ is a free-module over any 

104 system of parameters. Hence, if the set (771, . . . , rjt) is minimal for inclusion, K[x]*^ 

105 decomposes into a direct sum: 

t 

K[x]« = 077,.K[(?i,...,0„]. 

4 = 1 

106 This decomposition is called a Hironaka decomposition of the invariant ring. The 

107 9i are called primary invariants, and the rji secondary invariants (in algebraic com- 

108 binatorics literature, the 0i are some times called quasi-generators and the r]i sepa- 

109 rators |GS84j ). It should be emphasized that primary and secondary invariants are 

110 not uniquely determined, and that being a primary or secondary invariant is not an 

111 intrinsic property of an invariant p, but rather express the role of p in a particular 

112 generating set. 

113 The primary and secondary invariants together form a generating set, usually 

114 non minimal. From the degrees (di, . . . , d„) of the primary invariants (^i, . . . , On) 

115 and the Hilbert series we can compute the number t and the degrees (d'^, . . . ,d[) of 
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116 the secondary invariants (771, . . . ,rit) by the formula: 

z'^'^+--- + z< = (l-z'^i)---(l-z'*")iJ(K[x]<=,z). 

117 We denote this polynomial by 5(K[x]'^, z). Assuming di < • • • < rfn and d'l < • • • < 

118 dj, it can be proved that: 

d't^di + --- + d^-n- fi, /3(K[x]<=) < max(d„,4), 

|G| 

119 where /x is the smallest degree of a polynomial p such that a ■ p ~ det (cr)p for all 

120 (7 G G |Sta79[ Proposition 3.8]. 

For example, if G is the symmetric group 6„, the n elementary symmetric poly- 
nomials (or the n first symmetric power sums) form a system of parameters, i = 1, 

= and 771 = 1. This is consistent with the fundamental theorem of symmetric 
polynomials. More generally, if G is a permutation group, the elementary symmet- 
ric polynomials still form a system of parameters: K[x]'^ is a free module over the 
algebra Sym(x) = K[x]®" of symmetric polynomials. It follows that: 

121 For a review of algorithms to compute primary invariants with minimal degrees, 

122 see |DK02] . They use Grobner bases, exploiting the property that a set 0i, . . . , 8„ 

123 of n homogeneous invariants forms a system of parameters if and only if x = 

124 is the single solution of the system of equations 8i(x) — ■ ■ ■ = 0„(x) = (see 

125 e.g. pK02~t Proposition 3.3.1]). 

126 We focus here on the second step: we assume that primary invariants 0i, . . . , 8„ 

127 are given as input, and want to compute secondary invariants. This is usually 

128 achieved by using the following proposition to reduce the problem to linear algebra. 

129 Proposition 2.1. Let 0i, . . . , 0„ be primary invariants and S :~ (771, . . . , rjt) be a 

130 family of homogeneous invariants with the appropriate degrees. Then, the following 

131 are equivalent: 

132 (i) S is a family of secondary invariants; 

133 (ii) S is a basis of the gwotzeTit IK[x]'^/(8i, . . . , 8„)][j[x]g; 

134 (Hi) S is free in the quotient K[x]/(8i, . . . , 8„)k[x] • 

135 The central problem is how to compute efficiently inside one of the quotients 

136 K[x]/(ei, . . . , 8„)k[x] or IK[x]'^/(ei, . . . ,e„)K[x]G. Most algorithms rely on {Hi) 

137 using normal form reductions w.r.t. the Grobner basis for 81, ... , 8„ which was 

138 calculated in the first step to prove that they form a system of parameters. The 

139 drawback is that Grobner basis and normal form calculations do not preserve sym- 

140 metries; hence they cannot be used to confine the calculations into a small subspace 

141 of K[x]/(8i, . . . , 8„)k[x]- Besides, even the Grobner basis calculation itself can be 

142 intractable for moderate size input {n = 8) in part due to the large multiplicity 

143 [di - ■ ■ dn) of the unique root x ~ of this system. 

144 An other approach is to use (ii). Then, in many cases, one can make use of the 

145 symmetries to get a compact representation of invariant polynomials. For example, 

146 if G is a permutation group, an invariant can be represented as a linear combination 

147 of orbitsums instead of a linear combination of monomials, saving a factor of up to 

148 |G| (see e.g. |Thi01) ). Furthermore, one can use SAGBI-Grobner bases (an analogue 

149 of Grobner basis for ideals in subalgebras of polynomial rings) to compute in the 
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150 quotient (see |ThiOH IFR09] ). However SAGBI and SAGBI-Grobner basis tend to 

151 be large (in fact, they are seldom finite, see |TT04j ) . even when truncated. 

152 In both cases, it is hard to derive a meaningful bound on the complexity of 

153 the algorithm, by lack of control on the behavior of the (SAGBI)-Grobner ba- 

154 sis calculation. In the following section, we propose to calculate in the quotient 

155 K[x]'^/ (01, . . . , 0n)K[x]<^ using instead evaluation techniques. 

156 3. QUOTIENTING BY EVALUATION 

157 Recall that, in the good cases, an efficient mean to compute modulo an ideal is 

158 to use evaluation on its roots. 

159 Proposition 3.1. Let P be a system of polynomials in IK[x] admitting a finite 

160 set Pi , . . . , p,. of multiplicity-free roots, and let I be the dimension ideal they 

161 generate. Endow further W with the pointwise (Hadamard) product. Then, the 

162 evaluation map: 

$ : K[x] — > W 

P ' > {p{Pl),---,PiPr)) 

163 induces an isomorphism of algebra from K[x]//. In particular, K[x]// is a semi- 

164 .simple basic algebra, a basis of which is given by the r idempotents {pi)i=i,...,r which 

165 satisfy pi{pj) = Sij; those idempotents can be constructed by multivariate Lagrange 

166 interpolation, or using the Buchberger-Mdller algorithm |MB82j . 

167 This proposition does not apply directly to the ideal (9i, . . . , 8„) because it has 

168 a single root with a very high multiplicity di . . . dn. The central idea of this paper is 

169 to blowup this single root by considering instead the ideal . . . , On-i, ©n — e)j 

170 where e is a non zero constant, and then to show that the grading can be used 

171 to transfer back the result to the original ideal, modulo minor complications. This 

172 approach is a priori general: assuming the field is large enough, the ideal 0i, . . . , 0„ 

173 can always be slightly perturbed to admit di • • • d„ multiplicity-free roots; those 

174 roots are obviously stable under the action of G, and can be grouped into orbits. 

175 Yet it can be non trivial to compute and describe those roots. 

176 For the sake of simplicity of exposition, we assume from now on that G is a per- 

177 mutation group, that Oi, . . . , &„ are the elementary symmetric functions ei, . . . , e„, 

178 and that e = (—1)"+^. Finally we assume that the ground field K contains the n-th 

179 roots of unity; this last assumption is reasonable as, roughly speaking, the invariant 

180 theory of a group depends only on the characteristic of K. With those assumptions, 

181 the roots Pj take a particularly nice and elementary form, and open connections 

182 with well know combinatorics. Yet we believe that this case covers a wide enough 

183 range of groups (and applications) to contain all germs of generality. In particular, 

184 the results presented here should apply mutatis mutandis to any subgroup G of a 

185 complex reflection group. 

186 Remark 3.2. Let p be a n-th primitive root of unity, and set p '.— (1, p, . . . , p"^^). 

187 Then, ei(p) = ...=; e„_i(p) = and Cnip) = £• 

188 Proof. Up to sign, ei(p) is the i-th coefficient of the polynomial 



n-l 

(X^-l) = \{{X-p^). 

j=0 



□ 
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189 For a G ©„, write .= a ■ p the permuted vector. It follows from the previous 

190 remark that the orbit {Pf,)ae&r, of P gives all the roots of the system 

ei(x) = • • • = e„_i(x) = e„(x) - e = . 

191 Let X be the ideal generated by ei, . . . , e„_i, e„ — e in K[x], that is the ideal of 

192 symmetric relations among the roots of the polynomial X" — 1; it is well known 

193 that the quotient K[x]/I is of dimension n\. We define the evaluation map $ : p e 

194 K[x] (p(p^))cr as in Proposition 3.1 to realize the isomorphism from K[x]/I to 

195 £: = K®". 

196 Obviously, the evaluation of an invariant polynomial p is constant along G-orbits. 

197 This simple remark is the key for confining the quotient computation into a small 

198 subspace of dimension n!/|G|, which is precisely the number of secondary invariants. 

199 Let S'^ be the subalgebra of the functions in £ which are constant along G-orbits. 

200 Obviously, is isomorphic to where L is any transversal of the right cosets 

201 in &n/G. Let I'^ be the ideal generated by (ei, . . . , e„_i, e„ — e) in K[x]'^; as the 

202 notation suggests, it is the subspace of invariant polynomials in X. 

203 Remark 3.3. The restriction of ^ onK.[x]^, given by: 

$ : K[x]G — > 

P ' > {p{Pa))<reL 

204 is surjective and induces an algebra isomorphism between K[x]'^/Z'^ and . 

205 Proof. For each evaluation point p^, cr Cz L, set 

Pp. - J2 Pp. ' 

TEcrG 



207 
208 



206 where pp^ is the Lagrange interpolator of Proposition 3.1 Then, their images 
{^{Pp ))a-eL a-re orthogonal idempotents and, by dimension count, form a basis of 



209 From now on, we call evaluation points the family {p^)^^^. 

210 We proceed by showing that the grading can be used to compute modulo the 

211 original ideal (ci, . . . , Cn 

) , modulo minor complications. 

212 Lemma 3.4. LetG be a subgroup of&n andK. be afield of characteristic contain- 

213 ing a primitive n-th root of unity. Let S be a set of secondary invariants w.r.t. the 

214 primary invariants ei, . . . , e„, and write (S')k for the vector space they span (equiv- 

215 alently, one could choose a graded supplementary of the graded ideal (ei, . . . , e„) in 

216 K[x]'^J. Write Sd for the secondary invariants of degree d. Then, 

for0^d<n: $(K[x]G) ^Sdh) , 

ford^n: <i>(K[x]G) ^ ms^h) ® <J>(K[x]G_ J . 

217 In particular, $ restricts to an isomorphism from (S')k to . 

218 Proof. For ease of notation, we write the Hironaka decomposition by grouping the 

219 secondary invariants by degree: 

t dmax 

^[^f = nM^i. . . . , e„] - (^<i)KK[ei, . . . , e„] , 

220 where dmax is the highest degree of a secondary invariant. Then, using that 

$(ei) = • • • = $(e„_i) = O^G and $(e„) = l^c , 
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221 we get that <I'(K[ei, . . . , e„]) = $(K[e„]) = K.I^g, and thus: 

= $(K[x]<=) = J2 H{Sd)KMK[ei, e„]) = ^ ^SdW) , 

where, by dimension count, the sum is direct. Using further that e„ is of degree n: 
$(K[x]f ) = mSdh) + mSd-nhen) + mSd-2n)Kel) + ■■■ 

-*((5'd)K)©*((S'd_„)K) ® $((S'<j_2«)k) ©••• 

222 The desired resuh fohows by induction. □ 



223 In practice, this lemma adds to Proposition 2.1 two new equivalent characteri- 

224 zations of secondary invariants: 

225 Theorem 3.5. Let G C ©„ he a ■permutation group, take ei,...,e„ as primary 

226 invariants, and let S ~ {rji, . . . ,rjt) be a family of homogeneous invariants with the 

227 appropriate degrees. Then, the following are equivalent: 

228 (i) S is a set of secondary invariants; 

229 (iv) ^{S) forms a basis of ; 

230 (v) The elements of ^{Sd) are linearly independent in 8'^ , modulo the subspace 

0<j<d, n I d-j 

231 Furthermore, when any, and therefore all of the above hold, the sum in (v) is a 

232 direct sum. 



233 Proof. Direct application of Lemma 3.4 together with recursion for the direct sum. 

234 □ 

235 Example 3.6. Let G = ,43 = ((1,2,3)) be the alternating group of order 3. In 

236 that case, p is the third root of unity j, and K = Q(j) = Q ©q Q.j ®q Q-J^- We 

237 are looking for nl/\G\ — 2 secondary invariants, whose degree are given by the 

238 numerator of the Hilbert series: 



if(K[x]^,z) = 



3\{l-z)^ {l-z3)J (l-z)(l-z2)(l-z3) 

Simultaneously, the Gn-orbit o/(l,p, p^) splits in two G-orbits. We can, for exam- 
ple, take as evaluation points the two G-orbit representatives Pq = (l,p,p^) and 



241 pf^2 1) — (Pi IjP ); '"^'^ evaluation morphism is given by: 

$ : K[x]'^ — > K2 



P ' > (P(P()),P(P(1.2)))) 

For example, <i>(l) = $(63) = (1,1); whereas ^(ei) — $(62) — 0. Let us evaluate 



the orbitsum of the monomial x\x2 = x(2, 1, 0), using Remark 4-4 



O(x(2.1'0))(p^j) ^j((2, 1,0), (0,1,2)) ^-((2,1,0), (1,2,0)) _^ ^-((2,1,0), (2,0,1)) ^g^- ^ 
0(x(2.1."))(p^^ 2^)^j((2, 1,0), (1,0,2)) _^ ^•((2,1,0,(0,2,1)) ^ ^-((2,1,0), (2,1,0)) ^g^-2 
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That is $(o(xfa;2)) = 3.(j, j^). It follows that: 
<I>(K[x]C^)=K.(l,l) 
cI>(K[x]f) = <i>(K[x]«) = {(0,0)} 

$(IK[x]^) K.(l, 1) ® K.(3, 3) = K.$(l) ® K.$(o(x?a;2)) . 

242 In particular, 1 and o{x\x2) are two secondary invariants, both over K orQ; 

K[x]'^ = Syni(x) ® Sym(x).o(a;?X2) . 

243 We consider now the two extreme cases. For G ~ ©„, there is a single eval- 

244 nation point and a single secondary invariant 1; and indeed, $(1) = (1) spans 

245 $(K[x]'-') = K. Take now G = {()} the trivial permutation group on n points. 

246 Then, the eva luati on points are the permutations of (1, j, j^, . . . , In that 



247 case, Theorem 3.5 states in particular that the matrix {j^"^''^)m,a, where m and a 

248 run respectively through the integer vectors below the staircase and through ©„, is 

249 non singular. 

250 4. An ALGORITHM FOR COMPUTING SECONDARY INVARIANTS BY EVALUATION 

251 Algorithm [T] is a straightforward adaptation of the standard algorithm to com- 

252 pute secondary invariants in order to use the evaluation morphism $ together with 

253 Theorem 13.51 

254 For the sake of the upcoming complexity analysis, we now detail how the required 

255 new invariants in each degree can be generated and evaluated in the case of a 

256 permutation group. 

257 It is well known that the ring K[x] is a free Sym(x)-module of rank nl. It admits 

258 several natural bases over Sym(x), including the Schubert polynomials, the descent 

259 monomials, and the monomials under the staircase. We focus on the later. Namely, 

260 encoding a monomial m = x" in K[x] by its exponent vector a = (ai,...,a„), 

261 m is under the staircase if cki ^ n — i for all 1 ^ i ^ n. Given a permutation 

262 group G C 6„, a monomial m is canonical if m is maximal in its G-orbit for the 

263 lexicographic order: cr(m) ^lox rn, Vu € G. The following lemma is a classical 

264 consequence of the Reynolds operator being a K[x] '^-module morphism. 

265 Lemma 4.1. Let M be a family of polynomials which spans K[x] as a Sym(x)- 

266 module. Then, the set of invariants {R{m) \ m G M} spans K[x]'^ as a Sym(x) 

267 module. 

268 In particular, taking for M the set of monomials under the staircase, one gets 

269 that the orhitsums of monomials which are simultaneously canonical and under the 

270 staircase generate K[x]'-' as a Sym{x) -module. One can further remove non zero 

271 integer partitions from this set. 

272 Proof. Let p e K[x]'^ be an invariant polynomial, and write it as p = X^mgA/ fm'm', 

273 where the /,„ are symmetric polynomials. Then, using that the Reynolds operator 

274 i? is a K[x]'^-module morphism, one gets as desired that: 

p = R{p) = i?( ^ /„,to) = f^nRim) . □ 

275 Remark 4.2. The canonical monomials under the staircase can be iterated through 

276 efficiently using orderly generation |Rea78[ IMcK98j and a strong generating sys- 

277 tem of the group G |Ser03j ; the complexity of this iteration can be safely bounded 
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Algorithm 1 Computing secondary invariants and irreducible secondary invariants 
of a permutation group G, w.r.t. the symmetric functions as primary invariants, 
and using the evaluation morphism $. 

We assume that the following have been precomputed from the Hilbert series: 

• Sd'- the number of secondary invariants of degree d 

(this is the coefficient of degree d of S'(K[x]'^, z)) 

• Cd'- the dimension of dim $(K[x]^) 

(this is Sd if d < n and Cd-n + Sd otherwise) 
At the end of each iteration of the main loop: 

• Sd is a set Sd of secondary invariants of degree d; 

• Id is a set of irreducible secondary invariants of degree d; 

• Ed models the vector space 'l>(K[x]^). 
Code, in pseudo-Python syntax: 

def Secondarylnvariants(G) : 

for d G {0, 1, 2, . . . , deg(S'(]K[x]^, 2))} : 
U = {} 
Sd = {} 
if d ^ n : 

Ejj = Ed— ri #Defect of direct sum of Theorem |3.5| 
else : 

Ed = {0} 

# Consider all products of secondary invariants of lower degree 
for {r), rj') £ Sfc x I; with k + I = d : 

if ^rjTi') i Ed : 
Sd = SdU{W} 
Ed = Ed e]K.$(W) 

# Complete with orbitsums of monomials under the staircase 
for m £ CanonicalMonomialsUnderStaircaseOfDegree(d) : 

if dim Ed ea : 

break #A11 secondary invariants were found 
r\ — OrbitSum(m) 
if $(??) ^ Ed : 

Id = Id U{77} 

Sd = SdU{r,} 

Ed =Ed eK.$(?7) 
return ({So, Si, . . . }, {Iq, Ii, . . . }) 



278 above by 0(71!), though in practice it is much better than that (see Figures^ and[7} 

279 and [Nicllj for details). 

280 Remark 4.3. Let x" be a monomial. Then, evaluating it on a point requires at 

281 most 0{n) arithmetic operations in Z. Assume indeed that has been precomputed 

282 in K and cached for all k in 0, . . . ,n — 1; then, one can use: 

283 where a is written, in the scalar product, as a permutation of {0, . . . , ?i — 1}. 

284 Remark 4.4. Currently, the evaluation of the orbitsum o(x") of a monomial on 

285 a point is carried out by evaluating each monomial in the orbit. This gives a 
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Figure 1. This figure plots, for n < 10 and for each transitive 
permutation group G C ©„, the average computation time per 
canonical integer vectors below the staircase. If one ignores the 
two top groups for each n (which correspond to 2l„ and 6„ re- 
spectively), the worst case complexity is seemingly roughly linear 
in the size of the result: 0{nC{G)). 



286 complexity of 0{n\G\) arithmetic operations in Z (for counting how many times 

287 each appears in the result) and 0{n) additions in K (for expressing the result in 

288 Kj. This can be roughly bounded by 0{\G\) arithmetic operations in K. This bounds 

289 the complexity of calculating $(o(x")) on all points by y^OdGI) = 0{n\). 

290 This worst case complexity gives only a very rough overestimate of the average 

291 complexity in our application. Indeed, in practice, most of the irreducible secondary 

292 invariants are of low degree; thus Algorithm [T] only need to evaluate orbitsums 

293 of monomials m of low degree; such monomials have many multiplicities in their 

294 exponent vector, and tend to have a large automorphism group, that is a small 

295 orbit. 

296 Furthermore, it is to be expected that such evaluations can be carried out much 

297 more efficiently by exploiting the inherent redundancy (d la Fast Fourrier Trans- 

298 form). In particular, one can use the strong generating set of G to apply a divide 

299 and conquer approach to the evaluation of an orbitsum on a point. The complex- 

300 ity analysis and benchmarking remains to be done to evaluate the practical gain. 

301 Finally, the evaluation of an orbitsum on many points is embarrassingly parallel 

302 (though fine grained), a property which we have not exploited yet. 

303 5. Complexity analysis 

304 For the sake of simplicity, all complexity results are expressed in terms of arith- 

305 metic operations in the ground field K = Q(p). This model is realistic, because, in 

306 practice, the growth of coefficients does not seem to become a bottleneck; a possi- 

307 ble explanation for this phenomenon might be that the natural coefficient growth 
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Figure 2. This figure plots, for n < 10 and for each transi- 
tive permutation group G C 6„, the number C"(G) := C(G) — 
catalan(n) + 1 of canonical integer vectors below the staircase for 
G which are not non zero partitions, versus the number n!/|G| of 
secondary invariants. The dotted lines suggest that, in practice, 
n\l\G\ < C'{G) < {n\/\G\r-'. 



308 would be compensated by the pointwise product which tends to preserve and in- 

309 crease sparseness. We also consider that one operation in K is equivalent to n 

310 operations in Q. This is a slight abuse; however dimQQ(p) = (j){n) > 0.2n for 

311 n < 10000 which is far beyond any practical value of n in our context. 

312 Theorem 5.1. Let G be a permutation group, and take the elementary symmet- 

313 ric functions as primary invariants. Then, the complexity of computing secondary 

314 invariants by evaluation using Algorithm^is bounded above by 0{nl^ + n!^/|Gp) 

315 arithmetic operations in K. 

316 Proof. To get this upper bound on the complexity, we broadly simplify the main 

317 steps of this algorithm to: 

318 (1) Group theoretic computations on G: strong generating set, conjugacy classes, 

319 etc; 

320 (2) Computation of the Hilbert series of K[x]'^; 

321 (3) Construction of canonical monomials under the staircase; 

322 (4) Computation by $ of the evaluation vectors of the orbitsums of those mono- 

323 mials; 

324 ( 5) Computation of products ^{ri)^{ri') of evaluation vectors of secondary in- 

325 variants; 
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326 ( 6) Row reduction of the evaluation vectors. 

327 The complexity of ([T]) is a small polynomial in n (see e.g. [Ser03j ) and is negligible 

328 in practice as well as in theory. ^ can be reduced to the addition of c polynomials of 

329 degree at most (2) , where c < \G\ < n\ is the number of conjugacy classes of G (the 

330 denominator of the Hilbert series is known; the mentioned polynomials contribute 

331 to its numerator, that is the generating series of the secondary invariants); it is 

332 negligible as well. Furthermore, by Remark 4.2 (3]) is not a bottleneck. 



333 Using Lemma 4.1 and Remark |4.4[ the complexity of Q is bounded above by 

334 0{nl^) (at most 0{nl) orbitsums to evaluate, for a cost of 0(n!) each). 

335 For a very crude upper bound for ([s]), we assume that the algorithm computes 

336 all products of evaluation vectors of two secondary invariants. This gives (n!/G)^ 

337 products in £^ which is in 0{nl/G)^. 

338 Finally, in (|6]), the cost of the row reduction of ©(n!) evaluation vectors in S'^ 

339 is of 0(n!3/|Gp). □ 

340 This complexity bound gives some indication that the symmetries are honestly 

341 taken care of by this algorithm. Consider indeed any algorithm computing sec- 

342 ondary invariants by linear algebra in K[a;]/ Sym(x)"'' (say using Grobner basis or 

343 orthogonal bases for the Schur-Schubert scalar product). Then the same estimation 

344 gives a complexity oi 0{nl'^ /\G\) (reducing n\ candidates to get n!/|G| linearly inde- 

345 pendent vectors in a vector space of dimension n\). Therefore, for G large enough, 

346 a gain of |G| is obtained. 

347 That being said, this is a very crude upper bound. For a fixed group G, one 

348 could use the Hilbert series to calculate explicitly a much better estimate: indeed 

349 the grading splits the linear algebra in many smaller problems and also greatly 

350 reduces the number of products to consider. However, it seems hard in general to 

351 get enough control on the Hilbert series, to derive complexity information solely in 

352 term of basic information on the group (n, |G|, ...). Also, in practice, there usually 

353 are only few irreducible invariants, and they are of small degrees. Thus only few of 

354 the canonical monomial need actually to be generated and evaluated. 

355 It is therefore essential to complement this complexity analysis with extensive 

356 benchmarks to confirm the practical gains. This is the topic of the next section. 



357 6. Implementation and benchmarks 

358 Algorithm [ij and many variants, have been implemented in the open source 

359 mathematical platform Sage |S^09] . The choice of the platform was motivated by 

360 the availability of most of the basic tools (group theory via GAP |GAP99) , cyclotomic 

361 fields, linear algebra, symmetric functions, etc), and the existence of a community 

362 to share with the open-source development of the remaining tools (e.g. Schubert 

363 polynomials or the orderly generation of canonical monomials) [SCcOSj . Thanks 

364 to the Cython compiler, it was also easy to write most of the code in a high level 

365 interpreted language (Python), and cherry pick just those critical sections that 

366 needed to be compiled (orderly generation, evaluation). The implementation is 

367 publicly available in alpha version via the Sage-Combinat patch server. It will 

368 eventually be integrated into the Sage library. 

369 We ran systematic benchmarks (see Figure [6] and [6| , comparing the results with 

370 the implementation of secondary invariants in Singular [GPS981 IKin07b| . Note 

371 that Singular's implementation deals with any finite group of matrices. Also, it 
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372 precomputes and uses its own primary invariants instead of the elementary sym- 

373 metric functions. Therefore, the comparison is not immediate: on the one hand, 

374 Singular has more work to do (finding the primary invariants); on the other hand, 

375 when the primary invariants are of smah degree, the size of the result can be much 

376 smaller. Thus, those benchmarks should eventually be complemented by: 

377 • Calculations of secondary invariants w.r.t. the elementary symmetric func- 

378 tions, using Grobner basis using Singular and Magma; 

379 • Calculations of secondary invariants using Singular and Magma; 

380 • Calculations of secondary invariants w.r.t. the elementary symmetric func- 

381 tions, using SAGBI-Grobner basis (for example by using MuPAD-Combinat |Thil 

382 IHT040 . 

383 A similar benchmark comparing Magma [CP96| and MuPAD-Combinat is presented 

384 in [Thiol 'j Figure 1] (up to a bias: the focus in MuPAD-Combinat is on a minimal 

385 generating set, but this is somewhat equivalent to irreducible secondary invariants). 

386 This benchmark can be roughly compared with that of Figure |6] by shifting by a 

387 speed factor of 10 to compensate for the hardware improvements since 2001. Related 

388 benchmarks are available in |Kin07b| IKin07a| . 

389 We used the transitive permutation groups as test bed. A practical motivation 

390 is that there are not so many of them and they are easily available through the 

391 GAP database |Hul05j . At the same time, we claim that they provide a wide enough 

392 variety of permutation groups to be representative. In particular, the computation 

393 for non transitive permutation groups tend to be easier, since one can use primary 

394 invariants of much smaller degrees, namely the elementary symmetric functions in 

395 each orbit of variables. 

396 The benchmarks were run on the computation server sage . math . Washington . edr{^ 

397 which is equipped with 24 Intel(R) Xeon(R) CPU X7460 @2.66GHz cores and 

398 1 28 GB of RAM. We did not use parallelism, except for running up to four tests 

399 in parallel. The memory usage is fairly predictable, at least for the Sage imple- 

400 mentation, so we did not include it into the benchmarks. In practice, the worst 

401 calculation used 12 GB. Any calculation running over 24 hours was aborted. 

402 7. Further developments 

403 At this stage, the above sections validate the potential of the evaluation ap- 

404 proach. Yet much remains to be done, both in theory and practice, to design 

405 algorithms making an optimal use of this approach. The main bottleneck so far is 

406 the calculation of evaluations by $, and we conclude with a couple problems we are 

407 currently investigating in this direction. 

408 Problem 7.1. Construct invariants with nice properties under evaluation by $ 

409 (sparsity, ...). A promising starting point are Schubert polynomials |LS821 lLas03) . 

410 as they form a basis o/K[x] as SyTn{x) -module whose image under ^ is triangular. 

411 However, it is not clear whether this triangularity can be made somehow compatible 

412 with the coset distribution of G in 6„. 

413 Another approach would be to search for invariants admitting short Straight Line 

414 Programs. 



This server is part of the Sage cluster at the University of Washington at Seattle and is devoted 
to Sage development; it was financed by "National Science Foundation Grant No. DMS-0821725". 
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Figure 3. Comparative benchmark for the computation of sec- 
ondary invariants for aU transitive permutation groups for n < 10 
using Sage's evaluation implementation (+) and Singular's elim- 
ination implementation (o). The groups are sorted horizontally by 
increasing n and then by increasing cardinality. 



415 Note that a good solution to this problem, combined with the evaluation ap- 

416 proach of this paper, could possibly open the door for the solution of a long stand- 

417 ing problem, namely the explicit construction of secondary invariants; currently 

418 such a description is known only in the very simple case of products of symmetric 

419 groups |GS84) . Even just associating in some canonical way a secondary invariant 

420 to each coset in 6„/G seems elusive. 

421 From a practical point of view, the following would be needed. 

422 Problem 7.2. Find a good algorithm to compute $ on the above invariants. This 

423 is similar in spirit to finding an analogue of the Fast Fourier Transform w.r.t. the 

424 Fourier Transform. 
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Figure 4. Benchmark for the computation of secondary invariants 
for all transitive permutation groups for n < 10 (and for some 
below n < 14), using Sage's evaluation implementation. For each 
such group, n is written at position {k,t), where k = n\/\G\ is the 
number of secondary invariants and t is the computation time. In 
particular, the symmetric groups (3„ and the alternating groups 
2l„ are respectively above 1! and 2!. 



425 Theorem |3.5| further suggests that, using the grading, it could be sufficient to 

426 consider only a subset of the evaluation points. This is corroborated by computer ex- 

427 ploration; for example, for the cyclic group C7 of order 7, 110 evaluation points out 

428 of 720 were enough for constructing the secondary invariants. Possible approaches 

429 include lazy evaluation strategies, or explicit choices of evaluation points, or some 

430 combination of both. 

431 Problem 7.3. Get some theoretical control on which evaluation points are needed so 

432 that $ restricted on those points remains infective on some (resp. all) homogeneous 

433 component l&\yi\^ . 

434 Here again, Schubert polynomials are natural candidates, with the same difficulty 

435 as above. A step toward Problem |7.3| would be to solve the following. 

436 Problem 7.4. For G C 6„ a permutation group, and to start with for G the trivial 

437 permutation group, find a good description of the subspaces $(K[x]^). 
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442 
443 



Last but not least, one would want to generalize the evaluation approach to any 
matrix groups, following the line sketched in the introduction. The issue is whether 
one can get enough control on perturbations of the primary invariants so that: 

• The orbits of the simple roots are large, in order to benefit from the gain 
of taking a single evaluation point per orbit; 

• Only few of the primary invariants need to be perturbated, to best exploit 



the grading in the analogue of Theorem 3.5 
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